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This article considers the question of existence and uniqueness of the 
response of nonlinear time-varying RLC networks driven by independent 
voltage and current sources. It is proved that under certain conditions the 
response exists, is unique, and is defined by a set of ordinary differential 
equations satisfying some Lipschitz conditions. These conditions are of two 
types: (1) the network elements must have characteristics which satisfy 
suitable Lipschitz conditions and {2) the network must satisfy certain 
topological conditions. It shoidd be noted that elements with nonmonotonic 
characteristics are allowed and that the element characteristics need to be 
continuous but not differ entiable. 

I. INTRODUCTION 

This article considers the questions of existence and uniqueness of 
the response of nonlinear time- varying RLC networks. It is proved that, 
under conditions imposed on the network elements and the network 
topology the response exists, is unique, and is defined by a set of ordi- 
nary differential equations satisfying some Lipschitz conditions. Thus, 
from the conditions imposed on the network it follows that the response 
of a network of this class is continuous whenever the sources applied to 
the network arc continuous functions of time. In other words, for the 
class of networks under consideration, jump phenomena (of the type 
that occur in relaxation oscillators) are excluded. 1 

One motivation for studying this problem is the construction of non- 
linear net work models for physical devices and processes. The behavior 
of these models is often investigated by simulation studies performed on 
digital computers. It is clear that in order to get meaningful answers the 
existence and uniqueness of the model's response have to be assured. 
The simulation study requires the setting up of an appropriate set of 
differential equations and their integration. As networks of the class 

* On leave of absence from the Department of Electrical Engineering, Univer- 
sity of California, Berkeley, California. 
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considered here do not have jump phenomena, their equations can be 
integrated by some standard subroutines. 

This article may be viewed as an extension to the RLC case of the 
articles by R. J. Duffin 234 and G. Birkhoff and J. B. Diaz 5 which were 
devoted to nonlinear resistive networks. We make heavy use of topologi- 
cal considerations and had to extend the techniques developed for the 
linear case by many people 6 ' 7 - 8 P. R. Bryant in particular. 910 For further 
references see Ref. 16. 

In the next section, we classify the network elements and exhibit the 
basic assumptions which hold for the remainder of the article. Some 
simple nonlinear circuits are also considered. Section III presents some 
standard reductions of sources and the definition of determinateness. 
Section IV deals with one-element-kind networks; its theorems are 
generalizations of Dufnn's work and include some of his theorems as 
corollaries. The main result of the article is Theorem IV in Section V, 
which states the conditions under which a nonlinear RLC network is 
determinate. The conditions are of two types: (i) every characteristic 
has to satisfy suitable Lipschitz conditions and (it) the network has to 
satisfy certain topological conditions. It has to be noted that, first, 
elements with nonmonotonic characteristics are allowed and, second, 
that each characteristic has to be representable by a function which is 
continuous but not necessarily differentiable. Finally, in Section VI we 
introduce a symbolic notation which allows us to write the differential 
equations for the nonlinear case in a manner which resembles that of 
the linear case. 

II. ELEMENTS AND SIMPLE CIRCUITS 

2.1 Elements 

We assume that the reader has some familiarity with network theory, 
so that the basic concepts need not be defined. 1112 A network may be 
considered as a set of points, called nodes, and a set of connecting 
branches. Each branch represents a physical two-pole. We assume that 
the voltage drop across each two-pole and the current through each 
two-pole can be measured at any time. The sign conventions are shown 
in Fig. 1 : if, with respect to some arbitrary reference, the potential of A 
is larger (smaller) than the potential of B, then v is positive (negative) ; 
if the current actually flows in the direction of the arrow (opposite to the 
arrow) then i is positive (negative). Thus the product vi gives the power 
delivered by the outside world to the two-pole under consideration. 
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In most of the following, the branches consist of either a single source 
or a single element such as a resistor, an inductor or a capacitor. For each 
of these elements we shall adopt very broad definitions which we will 
narrow down in stating specific results. A two-pole is called a resistor 
if it is defined, for each t, by a set of ordered pairs (v, i), where v and i 
are finite numbers representing all the possible values, at time t, of the 
voltage and the current associated with the resistor. If the set of ordered 
pairs is independent of t, the resistor is said to be time-invariant. The set 
of (v, i) is called the characteristic of the resistor; for example, the charac- 
teristic of an ideal diode is given by 

{ (0,i): ^ i < «>}£/{ (vfi): - « < v ^ 0} . 

A resistor is called current-controlled if, for all time and all currents in 
the interval (—00,00), the voltage v(t) is a function* of the current i(t) 
and time I (we shall write u(t) = <R(i(0i0)> and the function <R(i,0 is 
a piecewise-continuous function f of / for each fixed number i. A voltage- 



v 
Fig. 1 — Sign conventions for two-pole. 

controlled resistor is defined in the dual manner. For example, a voltage 
source is a current-controlled resistor and a current source is a voltage- 
controlled resistor. If a resistor is current-controlled and time-invariant 
then the characteristic can be represented by a function v = 61 (i). A 
resistor is called a one-to-one resistor if, for each t, the voltage is related 
to the current by a one-to-one mapping from (— =0 , 00 ) onto (— 00 , co ) 
which may depend on time. 

A two-pole is called an inductor if it is defined, for each t, by a set of 
ordered pairs (<p,i) which represent the instantaneous flux and current 
associated with the inductor. The voltage across the inductor is given by 
d = dtp/dt. The current-controlled inductor, the flux-controlled inductor 
and the one-to-one inductor are defined as in the case of resistors. In the 
first two cases, if the elements are time-invariant, we shall write <p = £ (*) 
and i = T(<p), respectively. 

* Unless specifically indicated, we follow modern usage: each function is Bingle- 
valued; i.e., to each element of its domain it associates one and only one element 
of its range. , . . 

f A vector-valued function of time is said to be piecewise-continuous if it is 
continuous in every finite interval except at a finite number of points where it is 
discontinuous. At these points the function has a finite limit on the left as well 
as on the right. 
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A two-pole is called a capacitor if it is defined, for each t, by a set of 
ordered pairs (g,v) which represent the instantaneous charge and voltage 
associated with the capacitor. The current through the capacitor is given 
by i = dq/dt. The charge-controlled capacitor, the voltage-controlled capaci- 
tor and the one-to-one capacitor are defined as in the case of resistors. In 
the first two cases, if the elements are time invariant, we shall write 
v = SD(g) and g = G(v), respectively. 

Throughout the article we consider only elements whose characteris- 
tics can be represented, at all times, by a function defined on the interval 
(— 00,00 ). For example, Fig. 2(a), (b) and (c) represents the charac- 
teristics at time t of three time- varying resistors; we consider only re- 
sistors of the type shown in Fig. 2 (a) and (b), since they are current- and 






Fig. 2 — Characteristics at time / of three time-varying resistors: (a) and (b) 
are current- and voltage-controlled, respectively, while (c) is neither current- 
nor voltage-controlled. 
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voltage-controlled; these characteristics can be represented by 
u(t) = di(i(t),t), and i(t) = 9(f(0,0, 

respectively. The characteristics of Fig. 2(c) cannot be represented in 
this way, and resistors of this type will not be considered. 

Throughout the paper, whenever time-varying network elements are 
considered, it is assumed that the functions (R(-,0> S('.0, £(-,0, 
r (•,/), D (•,/), 6('i0 are piecewise-continuous functions of t for all 
fixed values of their first argument. 

In addition to resistors, capacitors and inductors, our networks in- 
clude voltage and current sources. Throughout this article we shall as- 
sume that the voltages of the voltage sources and currents of current 
sources are regulated functions of time.* For convenience we shall say 
that an element is continuous and monotonically increasing, when we 
mean that its characteristic is represented by a continuous monotonically 
increasing function which is defined on (— qo , oo ). 

It is convenient to refer to functions like (R(-,0 and 9D(*,i)i which 
represent the characteristics of some elements, as the characteristics of 
the elements. This slight misuse of the concept of a function and a rela- 
tion will be used only when there is no danger of confusion between the 
two. 

2.2 Two-Poles and Simple Connections of Two-Poles 

A two-pole is called voltage-controlled [current-controlled] if, for any 
initial time / and for any initial state, the voltage ?•(•) [the current i(-)] 
from /n on across its terminals determines uniquely the current i(-) 
through [the voltage »(■) across] the two-pole for / ^ / . 

A two-pole is said to be one-to-one if (a) it is both current-controlled 
and voltage-controlled and (b) it satisfies the following condition: for 
any initial state s , any initial time / , and any input current i(-), let 
/(sd ,i) be the voltage appearing at the terminals; for any initial state 
so , any initial time h , and any input voltage v(-), let g(s ,v) be the 
current — then it is required that 

g(so,f(so,i)) = i 

for all initial states .s and all input currents i{- ). 

An immediate consequence of these definitions is that any -parallel 
connection of a finite number of voltage-controlled two-poles is voltage-con- 
trolled. 



* A vector-valued function of time is said to be regulated when, for all t, it 
has a limit on the left as well as a limit on the right. 13 A step function and a rec- 
tangular wave are regulated functions. 
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Consider the case where there are only two two-poles in the parallel 
connection. Let them be characterized by the functions 

i k = F k (vMto)), (k = 1,2), 

where v is the voltage across the parallel connection, i k is the current 
through the /cth two-pole, s*(<o) is the state of the fcth two-pole at time 
to . The v*: and i* are real- valued functions defined on [/ , =° ). Kirchoff's 
current law implies that the current i through the parallel connection is 
given by 

Fi(v,8i (to)) + F 2 (v,s 2 (h)); 

hence, for fixed (si(/o),S2(fo))> i is a function of v. This argument obvi- 
ously extends, by induction, to the case where there are a finite number 
of two-poles. 

A dual argument would show that any series connection of a finite num- 
ber of current-controlled two-poles is current-controlled. 

A parallel connection of current-controlled two-poles is not necessarily 
current-controlled. Refer to Fig. 3, which shows the characteristics of two 
current-controlled resistors. The dashed line shows the characteristic of 
the parallel connection : depending on the operating point there may be 
three distinct values of the voltage for the same input current. Dually, 




CHARACTERISTIC OF THE 
/parallel COMBINATION 



CHARACTERISTIC OF 
SECOND RESISTOR 



CHARACTERISTIC OF 
FIRST RESISTOR 



Fig. 3 — Parallel connection of two current-controlled resistors. 
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SERIES CONNECTION 
OF g t AND g 2 



Fig. 4 — Series connection of two voltage-controlled resistors. 



a series connection of voltage-controlled two-poles is not necessarily voltage- 
controlled. 

To assume that each two-pole is one-to-one is not enough to cause both 
arbitrary parallel connections and arbitrary series connections to be 
one-to-one. Indeed, the well known characterization of continuous func- 
tions of bounded variation 14 implies that any voltage-controlled resistor 
characteristic, i(t) = (ii(v(t),t), that is continuous and of bounded variation 
in v can be obtained by connecting in 'parallel two one-to-one resistors whose 
characteristics are continuous and strictly monotonic. (One resistor is mono- 
tonically increasing and the other is monotonically decreasing.) A dual 
statement holds for current-controlled resistors. 

In fact, there are combined series and parallel connections of one- 
to-one two-poles that are neither voltage-controlled nor current-con- 
trolled. Refer to Fig. 4, which shows the series connection of Qj and Q 2 . 
Fig. 5 shows how a voltage-controlled characteristic such as Qi may be 
obtained by connecting in parallel two one-to-one resistors. Putting the 
two resistors of characteristic Qi and g 2 in series, we obtain (see Fig. 4) 
the characteristic g 3 , which is neither voltage-controlled nor current- 
controlled. 

A (possibly time-varying) flux-controlled inductor is a voltage-controlled 
two-pole and, dually, a (possibly time-varying) charge-controlled capacitor 
is a current-controlled two-pole. If the inductor is flux-controlled, the cur- 
rent i is a function of the flux ip:i(t) = T(<p(t),t). If v(- ) is the voltage 
applied to the inductor and <p is the flux through it at the initial time 
to , then by Lenz's law 
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*(f) = f v(t')dt' + 
J t 



<Po 



i(t) = r( / v(t')dt' + <po,t) for all t. 



■«) 



hence, 



2.3 Examples of One-To-One Two-Poles 

We present here a set of sufficient conditions under which some ele- 
mentary parallel or series connections of circuit elements constitute a 
two-pole which is either current-controlled, voltage-controlled or one-to- 
one. As the reader expects, quite specific assumptions will have to be 



CHARACTERISTIC OF 
FIRST ONE-TO-ONE 
RESISTOR 




CHARACTERISTIC OF 
SECOND ONE-TO-ONE 
RESISTOR 



Fig. 5— Parallel connection of two one-to-one resistors. 

made on the characteristics of the elements in order for the circuit to be 
a one-to-one two-pole. 

The elements that we are going to consider are capacitors, resistors 
and inductors. Let us rank order these elements together with voltage 
sources and current sources in the following way: E,C,R,L,J. We shall 
say that a resistor is higher in rank than an inductor or a current source 
but lower in rank than a capacitor or a voltage source. 

Until the end of this section, to simplify the discussion and without 
loss of generality, elements are assumed to be time-invariant. 

Theorem: Consider the following circuits: the parallel RC, the parallel 
RL, the parallel LC and the parallel RLC circuit. 

(A ) If (a) the highest-ranked element is current-controlled (charge- 
controlled in the case of the capacitor ) , 

(b) all other elements are voltage-controlled (flux-controlled in the case 
of the inductor), and 
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(c) the characteristics of all elements satisfy a Lipschitz condition ac- 
cording to Table I, 

then the parallel circuit is current-controlled. 
(B) If, in addition (d) the highest ranked element is one-to-one, 

then each parallel circuit is one-to-one. 

Proof: We shall consider only the RLC circuit, since the proofs of the 
simpler cases follow in a similar way. 

First let us prove that (a), (/;), and (c) imply that the circuit is cur- 
rent-controlled. Let i s be the source current. Then with the usual nota- 
tion 

*• = i + S(«) + r(*>) 

or, equivalently, using the fact that the capacitor is charge-controlled 

<i = -s(»fa)) - r(«0 +*.. (i) 

W = »(?) (2) 

By assumption (c) 33, T and £ satisfy Lipschitz conditions. Since the 
composite of two Lipschitzian functions is Lipschitzian, S(D(-)) is 
also Lipschitzian; therefore the system (1), (2) has a unique solution 
for each initial state and each current source. In this case the state is 
(q,<f). Thus the RLC is current-controlled. 

Second we prove that (a), (b), (c), (d) imply that the RLC circuit 
is one-to-one. It is immediate that these assumptions imply that the 
RLC circuit is voltage-controlled. It remains to show that it is one-to-one. 

Call <7i(-)> tpi(-), and i'i(-) the charge, flux and voltage resulting 
from the initial state (q n , <f>u) at time to and the input current i, . The 
functions qi ( ■ ) and fi(-) are the corresponding solutions of (1 ) and (2) ; 
i>i (t) = <p\(t) = £)(<7i (/.)). We have to show that, starting from the same 
initial state (f/ (l , vm) at time to , the input current resulting from the ap- 
plied voltage v, is precisely i K . 

Let q 2 , <?■> and i, be the resulting charge, flux and input current. It is 
immediate that v,(t) = ip-<(t) = £>(&(0)> Since <p 2 (to) = po> we have 

Table I 



Circuit 


Hinhest-Ranked Element 


Characteristics That Satisfy 
Lipschitz Conditions 


RL 
EC 
LC 
RLC 


R 
C 
C 
C 


(R(t'), rt» 

£>(9), rw 
©(?), 9(»), r(*) 
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<P2 = <fii • Since the capacitor is one-to-one, q 2 is uniquely denned by the 
relation above in terms of v\ , hence q$ = Qi . Finally, by Kirchhoff' s 
current law 

it = ?2 + 9( y z) + T(<p 2 ) 

= jfr + g(»(fc)) + r (*>,). 

The last expression is precisely i 8 by (1). Therefore i% = i, . This con- 
cludes the proof that the parallel RLC circuit is a one-to-one two-pole. 
The dual case is covered by the following 

Theorem: Consider the following circuits: the series RL, the series RC, 
the series LC and the series RLC circuit. 

(A) If 

(a) the lowest -ranked element is voltage-controlled (flux-controlled for 
inductor), 

(b) all other elements are current-controlled (charge-controlled for 
capacitors), and 

(c) the characteristics of all elements satisfy Lipschitz conditions ac- 
cording to Table II, 

then the series circuit is voltage-controlled. 

(B) If in addition (d) the lowest-ranked element is one-to-one, 
then each series circuit is one-to-one. 

The proof is similar to that of the previous theorem and is therefore 
omitted. 

III. REDUCTION OF THE NETWORK 

Throughout the article we consider networks consisting of nonlinear 
time-varying resistors, capacitors, inductors (without mutual induct- 
ance) and independent sources. We shall label by 91 the network under 
consideration. Usually, we consider each element and each source as 
constituting a branch of 31. We denote by E,C,R,L,.I the set of branches 
of 31 which are voltage sources, capacitors, resistors, inductors and cut- 
Table II 



Circuit 


Lowest-Ranked Element 


Characteristics That Satisfy 
Lipschitz Conditions 


RL 
RC 
LC 
RLC 


L 
R 
L 
L 


cR(t), rw 
£>(?), 9(») 
2) (q), W 

J)(?),(R(«),rw 
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rent sources, respectively. In our discussion, certain networks derived 
from 91 will play an important role. In order to refer to them conveni- 
ently, let us define the following notations: Let "A" be a subset of the 
set of branches of 91. Let us define 9 

9l x to be the network derived from 91 by removing all branches ex- 
cept the ones which are members of A, 

91(a) to be the network derived from 91 by replacing branches of set A 
by a short circuit, and 

9l (/( ). to be the network derived from 91 by removing the branches of 
set A. 

We shall use these notations as well as combinations of them. For ex- 
ample, 9l (B )c is the network derived from 91 by first replacing branches 
of set E (the voltage sources) by short circuits and then removing all 
elements which do not belong to set C. Similarly, 9l(*)(/>* is the network 
derived from 91 by shorting all voltage sources and removing all current 
sources. 

»S'*9l is defined to be the network derived from 91 by separating it into 
the maximum number of separable subnetworks. 

Throughout the article we assume that, first, for any cut set of current 
sources only, the source currents satisfy Kirchoff's current law, and, 
second, for any loop of voltage sources only, the source voltages satisfy 
Kirchoff's voltage law. 

Without loss of generality we consider networks that are connected 
and nonseparable. This assumption does not exclude the possibility 
that 9l(« )( y). be both unconnected and/or separable. In the following we 
shall prove that without a loss of generality we can restrict the discussion 
to a network 91 such that 9l(*)(j>« is both connected and nonseparable. 
The proof consists of an algorithm which changes the configuration and 
reduces 91 into a network 9l' which has the following properties: 

(0 9l' consists of connected subgraphs, 9l/, such that for each one of 
them, 9l'. (B)(J) . is connected and nonseparable. 

(ii) For all branches of 91 and 9l' which are not sources, any set of 
branch currents and voltages is a solution of 91 if and only if it is also a 
solution of 9l' (when the latter is driven by the corresponding sources). 

(ttt) Current sources of 9l' are linearly related to the current sources 
of 91. The same is true for voltage sources. 

The step-by-step reduction of the network 91 to 91 is done as follows. 

(1 ) From each loop which consists of voltage sources only, remove one 
voltage source. 

(2) In each cut set which consists of current sources only, replace 
one of the current sources by a short circuit. 
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The resulting network is connected; it has a tree which includes all the 
voltage sources as tree branches and all the current sources as links. 

(3) Each current source J whose fundamental loop includes more 
that one tree branch is removed from the network and is replaced by a 
set of current sources identical to J, each one placed in parallel with a 
tree branch of the fundamental loop. 

(4) All current sources that are in parallel with voltage sources are 
removed. 

(5) Any parallel connection of current sources is replaced by one 
equivalent current source. 

(6) Consider each fundamental cut set defined by a voltage source. 
For each one of them insert in every link a voltage source equal to that 
which is in the tree branch and, finally, short circuit the tree branch 
voltage source. 

(7) In each link, replace any series connection of voltage sources by 
one equivalent voltage source. 

(8) Separate the network into the maximum number of connected, 
nonseparable subgraphs. 

The resulting network is called 91 . Property (it) follows from the fact 
that all the steps of the above algorithm do not change the source con- 
tribution to any of the fundamental loop equations or the cut set equa- 
tions. Property (i) follows from the fact that all current sources are links 
of 91 and all voltage sources are in a link. Property (Hi) follows from 
steps (5) and (7). Finally, observe that jS*9l(j)(.r)* is identical with 
91 (*)(/)•. 

It is well known that the state of the network is completely determined 
by all the voltages, fluxes, charges, and currents in the branches of the 
network. In the case of linear networks it is well known that certain 
proper subsets of these variables may be chosen as the state. For special 
classes of nonlinear RLC networks similar subsets will be indicated in 
the sequel. 

We call a solution of an RLC network any set of voltages and currents 
of resistors, charges and voltages of capacitors, fluxes and currents of 
inductors which satisfy the Kirchhoff's laws and the branch characteris- 
tics. A network 91 is said to be determinate if for any value of the initial 
state So , given at any initial time t , and for any value of the sources 
E(-)i J(*)> there exists one and only one solution for t ^ to on some 
nonvanishing interval [to , t a ) • 

In the following section we shall describe a broad class of nonlinear 
RLC networks which are determinate. 
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IV. ONE-ELEMENT-KIND NETWORK 

The purpose of this section is to establish a set of sufficient conditions 
under which a nonlinear (possibly time-varying) resistor network driven 
by a set of independent current sources and voltage sources has, for all 
possible inputs, one and only one set of branch voltages and branch 
currents that satisfy Kirchhoff's laws. Conditions under which the 
solution satisfies a Lipschitz condition with respect to the sources are 
also given. 

The analysis of nonlinear resistor networks is almost identical with 
that of nonlinear capacitor networks or nonlinear inductor networks. 
Since the nonlinear resistors are the most flexible elements, we shall 
develop our analysis in terms of resistor networks. 

Let us start by making three preliminary remarks: 

(i) Given a resistor network together with an arbitrary distribution 
of current sources, it is always legitimate to assume that there are no 
cut sets of current sources only. (Dually, that there are no loops of volt- 
age sources only.) 

(ii) Any voltage source in series with a resistor may always be ab- 
sorbed into a suitably redefined branch characteristic. Refer to Fig. 6, 
where V\ and Vi are the node voltages of nodes 1 and 2 referred to the 
same datum. Let the current through the resistor be given by its charac- 
teristic g (v,t ) ; since g (v,t ) = g(vi — v«_ — e,t ) and since e ( • ) is a known 
function of time, we may introduce a new branch characteristic gn ( • , • ) 
specified at each instant of time by 

<7i2 (wi - v 2 ,t) A g(vi — v-i — e(t),t). 

In other words, the voltage source e has been absorbed into the time de- 
pendence of g\2 . A similar reasoning applies to a current-controlled re- 
sistor in series with a voltage source. 

The dual case can be taken care of in the same manner: in this case, 
a current source which is in parallel with either a voltage-controlled or a 
current-controlled resistor can be absorbed into the branch. 

Thus, without loss of generality, a network of nonlinear resistors and 
sources can be thought of as a network of nonlinear time-varying resist- 



,« fWV: KS> ^ 



Fig. 6 — Voltage source in series with resistor. 
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ors with the understanding that the sources have been absorbed in the 
branch characteristics. 

(Hi) Thus when, as in Theorems I and II below, we consider a net- 
work of nonlinear time-varying resistors, we include the case of a net- 
work made up of time-varying resistors and of independent sources. 
There is no loss of generality in considering only connected networks, 
since it amounts to considering successively each separate part of an 
unconnected network. 

We turn now to the statement of the main theorems. 

Theorem I (Existence and Uniqueness): Consider a connected nonsepara- 
ble network 91 of nonlinear (possibly time-varying) resistors. In case the 
resistor joining node a to node /3 is voltage-controlled, its characteristic is 
defined by the function g a p(-,-) such that g a p(v a (t) — v p (t),t) is the cur- 
rent flowing through it at time t from node a to node /3; here v a and vp are 
the node-to-datum voltages of nodes a and /3. Similarly, if this resistor is 
current-controlled, its characteristic is defined by the function r tt/3 (-,- ) such 
that r al }(i a fi(t),t) is the voltage difference between node a and node /3 at 
time t; here i a p is the current through the resistor measured positively if it 
flows from a to /3. 

// 

(a) there exists a tree 3 such that all its tree branches are current-con- 
trolled and all its links are voltage-controlled, 

(b ) for all a,fi, all t and all x in ( — oo , oo ) 
g a p(x,t) = -gp a (-x,t) if (a,/3) is a link 
r a0 (x,t) = —r a0 (—x,t) if (a,/3) is a tree branch 

(c) for all links and all t, g a ft(-,t) is a monotonically (not necessarily 
strictly) increasing continuous function defined on (— «•,»), and for all 
tree branches and all t,r a p(-,t) is a monotonically (not necessarily strictly) 
increasing continuous function defined on (— °° , °° )• 

Then, 

for all current-sources i" connected between any pair of nodes and for 
all voltage sources e' connected in series with network branches there exists 
one and only one set of branch voltages and branch currents that satisfy the 
Krichhoff laws and the branch characteristics. 

The conclusion of Theorem I can also be stated as follows: any net- 
work 9l', formed from 91 by inserting any set of voltage sources in series 
with any branch and any set of current sources between any node pair, 
is determinate. 

Assumption (b) is a consequence of the physical meaning of the func- 
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tions g and r and of the sign conventions: from a physical point of view 
they do not restrict the nonlinear resistors in any way. The two corol- 
laries that follow are special cases of Theorem I. Corollary I is an ex- 
tension of Theorems 2 and 3 of Duffin, 3 and is implied by his 1948 pa- 
per. 4 Such an extension has been pointed out by I.W. Sandberg. 

Corollary 1: Consider a connected network of nonlinear voltage-con- 
trolled (possibly time- varying) resistors. 

// 

(a) for all branches and all t, g aP (-,t) is a, monotonically (not neces- 
sarily strictly) increasing, continuous function defined on (— »,<» ), and 

(b) there exists a tree 3 such that all its branches have g a p's which are, 
for all t, monotonically increasing one-to-one mappings of (—00,00) 
onto ( — oo , oo ) , 

then the conclusion of Theorem I holds. 

Proof: The conclusion follows directly from Theorem I since the tree 
branches have g a g's that are, for all /, monotonically increasing one-to- 
one mappings of ( — « , oc ) onto (— °° , °° ) ; hence the tree branches are 
also current-controlled resistors satisfying assumption (c) of Theorem 
I. 

Corollary 2: Consider a connected network of nonlinear, current-con- 
trolled (possibly time-varying) resistors. 

// 

(a) for all branches and all /, r a p(-,t) is a monotonically (not neces- 
sarily strictly) increasing, continuous function defined on (—00,00), and 

(b) there exists a tree 3 such that its links have rp a s which are, for 
all t, monotonically increasing one-to-one mappings of (—00,00) onto 
(— 00,00 ), then the conclusion of Theorem I holds. 

We consider now the extension of the Thevenin and Norton equiva- 
lent circuits to nonlinear resistive networks. If we pick an arbitrary node 
pair of such a network 91, we may regard these nodes as the terminals 
of a two-terminal network: we shall call the characteristic of this two- 
terminal network the input characteristic of 91 at these two nodes. Dually, 
if we pick a branch and insert two terminals in series with it, we obtain a 
two-terminal network: we shall call the characteristic of this two-ter- 
minal network the branch-input characteristic of 91. 

Theorem II {Thevenin and Norton equivalent circuits): Consider a net- 
work 91 satisfying the requirements of Theorem I together with the same 
kind of source distribution. 
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Then 

(a) the input characteristic of 91 at any node pair is that of a current- 
controlled resistor whose characteristic is a continuous, monotonically in- 
creasing function defined on (—00,00). This characteristic may be repre- 
sented by the Thevenin equivalent circuit of Fig. 7(a): a series combination 
of a voltage soiwce and a monotonically increasing current-controlled re- 
sistor whose characteristic passes through the origin. 

(b) The branch-input characteristic of 91 at any branch is that of a volt- 
age-controlled resistor whose characteristic is a continuous, monotonically 
increasing function defined on (— 00,00 ). This characteristic may be repre- 
sented by the Norton equivalent circuit of Fig. 7(b): a parallel combination 
of a current source and a monotonically increasing voltage-controlled resistor 
whose characteristic passes through the origin. 

Let us consider some special cases of Theorem II. 

Corollary 3: Consider a connected network of nonlinear (possibly time- 
varying) resistors satisfying assumptions (a), (b) and (c) of Theorem I. 

(a) If, in addition, the characteristics of the tree branches of 3 are 
strictly increasing, then the input characteristic at any node pair is 
that of a strictly increasing current-controlled resistor. If the character- 
istics of all tree branches of 3 are continuous, monotonically increasing, 
one-to-one mappings of ( — 00 , oc ) onto ( — =o , 00 ) then so is the input 
characteristic at any node pair. 

(b ) If the characteristics of the links of the tree 3 are strictly increas- 
ing, then the branch-input characteristic is that of a strictly increasing 
voltage-controlled resistor. If the characteristics of all links of 3 are 
continuous, monotonically increasing, one-to-one mappings of ( — 00 , 00 ) 
onto (— 00 , 00 ), then so is any branch-input characteristic. 

Proof of Theorems I and II: The proof of these two theorems is divided 
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Fig. 7 — (a) Thevenin equivalent circuit: series combination of voltage source 
and a monotonically increasing current-controlled resistor whose characteristic 
passes through the origin, (b) Norton equivalent circuit : parallel combination of 
current source and monotonically increasing voltage-controlled resistor whose 
characteristic passes through the origin. 
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into two parts: in part one, we show that if Theorem I holds for a /.--node 
network then Theorem II is true for a A-node network. In part two, we 
use this implication to prove Theorem I by induction. 

The statement of the theorem allows time-varying resistors (and 
hence includes independent sources); however, in order to have as sim- 
ple a notation as possible, we write down the proof as if all resistors were 
time-invariant. 

Part One: We show that, for any integer k ^ 2, if Theorem I holds for 
a A-node network then the input characteristic at any node pair is that 
of the Thevenin equivalent circuit specified in Theorem II (a). Let us 
connect the node pair under consideration to a current source i, (see 
Fig. 8); this current source is viewed as an additional link, since it is a 
voltage-controlled resistor. By assumption, to each i„ there is one and 
only one set of branch currents and voltages that satisfy Kirchhoff's 




+ 
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Fig. 8 — Node pair connected to current source. 

laws and the branch characteristics. Consider two distinct values of i„ , 
namely, i„ and i/. Let the corresponding branch variables be v,i and 
v',i'. For each current-controlled branch define a number T (which de- 
pends on i and i' ) by the relation 

v - v AAv = r(i) - r(i') A r-(i - i) = f-Ai. 

Since all the current-controlled branches are monotonically increasing, 
r ^ 0. (If Ai = 0, P may be taken to be any nonnegative number.) 
Similarly, we define a g for each voltage-controlled resistor; again g ^ 0. 
The set of Ay's and Ai's together with Av„ and Ai, may be considered as 
a set of branch voltages and branch currents together with the source 
voltage and source current of a linear resistive network which is obtained 
by replacing each current-controlled resistor by a linear resistor of re- 
sistance f, each voltage-controlled resistor by a linear resistor of con- 
ductance g and the current source by a current source Ai, . Since the 
Ay's and Ai's satisfy Kirchhoff's laws, Tellegen's theorem 12 holds, 

Av„ • Ai K = 22 AvAt 

where the sum is over all resistive branches. 
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Since all branches have monotonically increasing characteristics, this is 
a sum of nonnegative terms and Av s Ai s ^ 0. In other words, Ai, > 
implies that Av a ^ : that is, the Thevenin equivalent circuit has a cur- 
rent-controlled monotonically increasing characteristic. The continuity 
of the characteristic follows from the following considerations: irrespec- 
tive of the values of the f's and g's, the fact that assumption (c) of 
Theorem I requires them to be nonnegative implies that the current 
transfer ratio from the current source to any branch has a magnitude no 
larger than unity; 11 hence Ai s — » implies Ai — > for all branches. Since 
the tree branches have continuous characteristics, it follows that, for 
them, Ay — * and, by Kirchhoff's voltage law, the same holds for the 
links. Hence Ai, —* implies Av s —* 0, i.e., the current-controlled charac- 
teristic of the Thevenin equivalent circuit is continuous. The proof of 
part (b) of Theorem II follows exactly the dual of the above argument. 
Part Two: Let us prove Theorem I for a two-node network (see Fig. 
9). Let us plot on the (v,i) plane of Fig. 10 the characteristics of the cur- 
rent-controlled tree branch and that of the voltage-controlled link, tak- 
ing into account the sign conventions defined on Fig. 9. By assumption, 
the functions g and r are both continuous and have (— »,«) as do- 
mains; therefore their representative curves intersect at least at one 
point (v,i). We assert that it is the only one: indeed, suppose there were 
a second one, (v ,i ) ; then the nionotonicity of r and g imply, respectively 

(V - v) (i - i) ^ and (»' - v)(i - i) ^ 0. 

Hence 

(v - v) (t -0=0. 

Suppose v = v; then since g is a function 

i = g(—v) = g(-v') = i. 

Similarly, if i = i, the fact that r is a function implies v = v . Hence for 
all possible sources, there is one and only one set of branch voltages and 
currents that satisfies Kirchhoff's laws and the branch characteristics. 



-L=g(v)< v <r(l) 



Fig. 9 — Two-node network. 
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Fig. 10 — Characteristics of current-controlled tree branch and of voltage- 
controlled link as function of the tree branch current and voltage. 

Thus Theorem I is established for a two-node network. The next step in 
the proof of Theorem I is to show that if it is true for an n-node network 
it is true for an (n + 1 )-node network. Consider the n-node network 
shown in Fig. 11. We shall build out this network into an (n + l)-node 
network. 

Let us first connect the tree branch between node n and node (n + 1), 
i.e., a current-controlled resistor. (There is no loss of generality in as- 
suming that the numbering of the nodes is such that the branch (n,n + 1 ) 
is a tree branch.) It is obvious that, for this network, the existence and 
uniqueness of the solution holds for all sources. Consequently, from part 
one of the proof, the input characteristic at any two nodes of this par- 
ticular (n + l)-node network has the equivalent circuits specified by 
Theorem II. The next step is to add a link, say between node k and node 
(n + 1). Since the input characteristic at the node pair (k,n + 1) is as 
specified in Theorem II (a), the voltage and current in the link are 
uniquely determined by the reasoning given for the case n = 2, and 
consequently the distribution of voltages and currents in all branches of 
the network is uniquely determined. 



n-NODE 
NETWORK 



(n) 



(k) 




(n+i) 



(0 



Fig. 11 — iV-node network built out to (n + l)-node network. 
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The process of constructing the (n + 1 )-node network from the n-node 
network can be carried out step by step, adding a link at a time. Thus 
at the end of the process there is one and only one set of branch voltages 
and currents in the (n + l)-node network that satisfies Kirchhoff's 
laws and the branch characteristics. Q.E.D. 

For the purpose of solving the network differential equations of a gen- 
eral nonlinear RLC network it is important to know, for the resistive 
network case, under what conditions the function which maps the 
sources, (EJ), into the branch voltages and currents, (v,i), satisfies a 
Lipschitz condition. It is immediately clear that additional assumptions 
are required: consider Fig. 12, which shows the characteristic of a cur- 
rent-controlled resistor which fulfills the conditions of Theorem I. In 
the neighborhood of the operating point A, this resistor may appear to 
small signals either as an open circuit or a short circuit. Note that the 
same statement would apply if the resistor were voltage-controlled. It is 
obvious that under such conditions, the mapping (E,J) — * (v,i) will not 
satisfy a Lipschitz condition. As shown in the following theorem, only 
weak additional assumptions are required. 

Theorem III: Consider a connected network of nonlinear (possibly time- 
varying) resistors which satisfies conditions (a), (b) and (c) of Theorem 
I. If, in addition, the following Lipschitz conditions are satisfied: there is a 
real-valued function h (R,t ) , defined and positive for R > and all t, such 
that 

I g a p(x,t) - g a p{x,t) | ^ h(R,t) | x - x | 

for all links of 3, for all x, x in ( — R,R) and all t and 




Fig. 12 — Characteristic of current-controlled resistor fulfilling conditions of 
Theorem I; note unbounded slope at point A. 
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| r afl {xJL) - r a0 (x',t) \ ^ h(R,t) \x - x'\ 

for all branches of 3, for all x,x in ( — R,R) and all t, then the mapping 
which maps (E,J) into (v,i) satisfies a Lipschitz condition* 
Proof: Consider the effect of a change in the voltage sources E on the 
branch voltages v and branch currents i. E is the vector whose z'th com- 
ponent is the output voltage e, of the source located in the z'th branch. 
In the present (n + 1) node network there are at most n(n + l)/2 
branches, hence E has at most that many components. Suppose that the 
change from E to E -f- AE is obtained by changing e, to e, + Ae, suc- 
cessively with i — 1, 2, •• • . Call a,/3 the terminals of the first branch 
and call Ni the remainder of the network (see Fig. 13). Since the input 
characteristic of N\ is monotonically increasing, and since an increase of 
the voltage across the nonlinear resistance R increases the current 
through it or keeps it constant, the change in the input voltage of N\ , 
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Fig. 13 — Nonlinear resistor R and voltage source in one branch of (n + 1)- 
node network; A'i represents remainder of network. 



(i) 



Av a p due to the change of e t to e\ + Aei is such that | Aci | §£ | Av a p 
(The superscript 1 indicates that only the source voltage in the first 
branch has been changed.) Call At»*. (I) the corresponding change in the 
A-th branch voltage. We assert that 



AtV 



^ I At. a< 



^ I Ad 



For the particular change in the sources under consideration, we may 
define, as in the proof of Theorem I, for each tree branch a suitable t 
and for each link a suitable g. Observe now that A<?i and the A% ' may 
be interpreted as being the source voltage and the resulting branch volt- 
age of a linear network which has the same configuration as the given 
nonlinear network but in which each nonlinear resistor is replaced by 
f or g as required. By assumption (c) of Theorem I, all the i" s and g's 
are nonnegative, hence all the voltage transfer ratios | Aw* ! /Ae l | of 



* Incidentally, if the network 91 was derived from another network 91^ , by 
applying to 91 ,i the algorithm of Section III, then the mapping (E,J) — » (v,i) is 
one-to-one. 
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the linear network cannot 11 exceed 1 and the inequality asserted above 
follows. Thus, for all i's and fc's, 

I W° I ^ I Ae,-|. 

Let Avk be the change in the voltage across the Mb. branch when E be- 
comes E -f AE. Summing over i, using the triangle inequality, and de- 
fining the norm of a vector as the sum of the magnitude of its components, 
we get 

I A* I ^ ||AE||. (3) 

Since there are at most n(n + l)/2 branches, we get finally 

|| Av || =S Kn+ 1)/2]|| AE|| (4) 

where Av is the change in the branch voltages corresponding to the 
change of the voltage sources from E to E + AE. We next bound the 
change in the branch currents. Applying (3) to a link and using the 
Lipschitz condition we find that 

I Ai k | ^ h(R,t) || AE || (for all links) 

and since there are at most n (n — 1 )/2 links and the change in a tree 
branch current is equal to the change in the sum of currents of the links 
which belong to its fundamental cut set, 

I Ait | ^ h(R,l)[n(n - l)/2] || AE || (for all branches). 

Thus 

(I Ai|| ^ h{R,t)[n{n - l)/4] || AE ||. (5) 

The effect of a change in the current sources from J to J + AJ is ob- 
tained in a dual manner. Since the current transfer ratio may not exceed 

unity 11 we get 

I Ai, | £ || AJ II 
and 

||Ai|| :g [n(»+ 1)/2]||AJ||. (6) 

This implies 

I Avk | ^ h(R,t) || AJ || (for all tree branches) 

and, by Kirchhoff's voltage law, 

I Avii | ^ h(R,t)n || AJ || (for all branches). 
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Finally 

|| Av || ;g h(R,t)[n 2 (n 4- l)/2] || Aj ||. (7) 

Using the usual product topology 17 for both the product spaces of voltage 
sources and current sources on the one hand and branch voltages and 
branch currents on the other, and invoking (4) to (7), we conclude that 
the mapping (EJ) — * (v,i) is Lipschitz. 

V. NONLINEAR RLC NETWORKS 

The previous section required all elements of the network to be of the 
same kind and to have a monotonically increasing characteristic. In this 
section both requirements are removed. In addition to independent 
sources, the network consists of nonlinear (possibly time-varying) 
resistors, capacitors and inductors and some of the elements are allowed 
to have characteristics with negative slope. 

As a first step let us make one remark. Theorems I, II, and III would 
still hold if all resistors were montonically decreasing instead of mono- 
tonically increasing. In the more complicated situation considered here 
the same possible choice exists. For example, separable subnetworks of 
9l<s)<7 which contain more than one capacitor could just as well contain 
monotonically decreasing capacitors. For simplicity, we shall assume 
that all monotonic elements are increasing. 

In order to state the following theorem we need two definitions. A 
network (or subnetwork) is called a self -loop if it consists of a single 
branch whose end-points are identified: it consists of one branch and 
one node. A network (or subnetwork) is called an open branch if it con- 
sists of a single branch whose end-points are not identified: it consists 
of one branch and two nodes. 

Theorem IV: Let 91 be a network of independent sources and nonlinear 
(possibly time-varying) resistors, capacitors and inductors (without mutual 
inductance) such that: capacitors of 91 are either charge-controlled or mono- 
tonically increasing voltage-controlled; resistors are either voltage-controlled 
or current-controlled; inductors are either flux-controlled or monotonically 
increasing current-controlled. It is further assumed that 91 and 9l (S )(j)» are 
nonseparable and connected. The network 91 is determinate if: 

(1) The capacitor network £*9l( B) c satisfies the following requirements: 

(a) Open branches of S*VL {E )c are charge-controlled and contain all 
charge-controlled capacitors which are not monotonically increasing. 

(b) Each subnetwork of .S*9l( B) c which contains more than one element 
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has a tree with monotonically increasing charge-controlled tree branches and. 
monotonically increasing voltage-controlled links. 
(c) Self -loops of S*3L(e)c are voltage-controlled. 

(2) The resistive network *S*91(.ec)r satisfies the following requirements: 

(a) Open branches are current-controlled and contain all current-con- 
trolled resistors which are not monotonically increasing. 

(b) Each subnetwork which contains more than one element has a tree 
with monotonically increasing current-controlled tree branches and mono- 
tonically increasing voltage-controlled links. 

(c ) Self -loops are voltage-controlled and contain all voltage-controlled re- 
sistors which are not monotonically increasing. 

(3) The inductive network S*?fl(scR)L satisfies the following require- 
ments: 

(a) Open branches are current-controlled. 

(b) Each subnetwork which contains more than one element has a tree 
with monotonically increasing current-controlled tree branches and mono- 
tonically increasing flux-controlled branches. 

(c) Self -loops are flux-controlled and contain all flux-controlled in- 
ductors that are not monotonically increasing. 

(4) In any finite interval, and for all time, the characteristics of the net- 
work resistors, capacitors and inductors satisfy a Lipschitz condition with 
respect to the following variables: 

tree branches: capacitors, with respect to q 

resistors and inductors, with respect to i 
links: capacitors and resistors, with respect to v 
inductors, with respect to <p. 

Remarks: Note that nonmonotonic voltage-controlled capacitors and 
current-controlled inductors were excluded from the discussion. Such 
capacitors and inductors may be included in the discussion provided they 
fall into the following trivial cases: each nonmonotonic voltage-con- 
trolled capacitor is in parallel with a voltage source and each nonmono- 
tonic current-controlled inductor is in series with a current source. In 
such cases, 9l(*) ( j). is separable unless 91 contains one element only. 

The above conditions insure the existence and uniqueness 18 of the 
solution on some nonvanishing interval [to , l a ), where t a > t Q . The 
length of this interval cannot be specified without further assumptions 
on the Lipschitz constants h(R,t). This is the well known problem of 
finite escape time. In particular, if for all branch characteristics the same 
Lipschitz constant can be used and holds over the whole domain of the 
characteristic, then the solution exists and is unique on [to , °o ) for all 
regulated E's and J's. 
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Proof of Theorem IV: Let us denote the voltages and charges of the ea- 
pacitive branches by (e c ,q r ). Similarly, denote the voltages and cur- 
rents of the resistive branches by (e r , i r ) and fluxes and currents of in- 
ductive branches by (*q l ,i L ). Voltage sources and current sources will 
denoted as usual by E and J. 

We assert that conditions (2) and (4) of Theorem IV imply, first, 
that the currents and the voltages of the resistive branches at time t, 
(e r (0,ir(0) are uniquely determined by the values, at the same time t, 
of the capacitor voltages, the voltage sources, the inductor currents and 
the current sources, (e c (t),E(t),i L (t)J(t)); and, second, that the map- 
pings 

e r (0 = f er (e E (t),E(t),i L (t)J(t),t) (8) 

ir(0 = Me t (0,E(0,UOJ(0,0 (9) 

satisfy Lipschitz conditions. 

Given any set of capacitor voltages e c and inductor currents i L such 
that Kirchoff's voltage law is satisfied in each loop formed by capacitors 
and voltage sources, and such that Kirchoff's current law is satisfied in 
each cut set formed by inductors and current sources, let us replace each 
capacitor with a voltage source whose voltage is equal to the voltage of 
the replaced capacitor and replace each inductor with a current source 
whose current is equal to the current of the replaced inductor. The net- 
work consists now of resistors, current sources and voltage sources only. 

Let us use the algorithm of Section III to change the configuration of 
the sources and to separate the network into its separable parts. Let us 
denote the resulting network by Tfl" and its sources by (E ff ,J*). 

The network 31* has three sets of subnetworks: (a) connected non- 
separable subnetworks which contain sources and two or more resistive 
branches, (b) subnetworks containing one resistive branch in parallel 
with a voltage source, and (c) subnetworks containing one resistive 
branch in parallel with a current source. 

Consider the first set of subnetworks. Denote the branch voltages and 
branch currents of these subnetworks by (e r , i r )i , and their sources 
by (E/'Ji"). From conditions (2) and (4) of Theorem IV it follows that 
each subnetwork contains a tree whose resistive branches are mono- 
tonically increasing current-controlled resistors and whose links are 
monotonically increasing voltage-controlled resistors, and that all ele- 
ments satisfy Lipschitz conditions. Therefore, from Theorems I and III 
of Section IV it follows that (e r , i r )i are uniquely denfied by (Ei",Ji R ) 
and that the mapping (e r (0,i<-(0) = f r (Ei" (t), J i" (t),t) satisfies Lip- 
schitz conditions. 
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The second set of networks corresponds to self -loops of S* l 3l iBC )R . In 
each subnetwork, the resistor is voltage-controlled, and hence the cur- 
rent through it is uniquely determined in terms of the voltage source. 
Since the characteristics satisfy a Lipschitz condition, the mapping from 
space E 2 B (the voltage sources) to (e r , i r )2 (the branch voltages and 
currents in the subnetworks of the second set ) satisfies Lipschitz condi- 
tions. 

Subnetworks of the third set correspond to open branches of &*9l( BC )« . 
As each subnetwork contains only one resistor and a current source, the 
requirement that the branch be current-controlled is enough to insure 
uniqueness of (e r , i r )3 , the branch voltages and currents of this set, in 
terms of the corresponding sources J 3 *. From condition (4) it follows 
that the mapping from J 3 R to (e r , i r )3 satisfies the Lipschitz condition. 

The voltages of sources E fi are linear combinations of the voltages E 
and e c , and the currents J R are linear combinations of J and i L (see 
Section IV). From this linearity property and from the properties of the 
above relations between the voltages and currents of the resistive 
branches and (E K J R ) it follows that (e r (04r(0) are uniquely defined 
by (e c (0,E(0,u(0 J(0) and that the mappings in (8) and (9) satisfy 
Lipschitz conditions. 

Let us now consider the capacitors of the network 91. Given any set of 
resistor currents i r and inductor currents i L such that Kirchoff 's current 
law is satisfied in each cut set formed by resistors, inductors and current 
sources, let us replace the inductive and resistive branches of 91 by 
current sources with currents equal to the corresponding currents i L , 
i r . The network consists now of sources and capacitors only. Let us use 
the algorithm of Section III to change the configuration of the sources 
and separate the network into its separable parts. The resulting network 
is denoted by 91 c and its sources by E C ,J C . We are going to establish an 
analogy between 91 c and its sources by E C ,J C . We are going to establish 
analogy between 91 c and 91 R and use the result just proved for 91" to 
deduce a similar result for 91 c . 

9l c consists of the three sets of subnetworks which were described in 
connections with 91* . Consider the second set of subnetworks of 9l c , 
which consists of single capacitors in parallel with a voltage source. 
Except for the trivial case where 91 consists only of a single capacitor in 
parallel with a voltage source, this set is empty, for otherwise 9l(B)(y>» 
would be separable. Condition (1 ) implies that each subnetwork of the 
first set has a tree, say r c , whose tree branches are monotonically increas- 
ing and charge-controlled, and whose links are monotonically increasing 
and voltage-controlled. For each subnetwork of this set, with each funda- 
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mental cut set of r c denned by a capacitive tree branch, we assign a 
variable 5, equal to the sum of the charges on all the capacitors of that cut 
set. For each subnetwork of the third set, we assign a q t equal to the 
charge on the capacitor, q will denote the vector whose components are 
the qi's. 

The analogy between 91 c and 91* is established in four steps: 
(i) For 91*, the Kirchhoff current law applied to the ith cut set as- 
sociated with a resistive tree branch reads 

./,* = £*" 

k 

where i ki is the current in the fcth branch of the ith cut set. The i ki are 
components of i r . For 9l c , we have by definition of q { , 

k 

where q ki is the charge in the kth branch of the ith cut set. The q k i are 
components of q. 

(m) For both 91* and 9l c , the Kirchhoff voltage law holds. 

(in) Condition (1) imposes requirements on the topology and ele- 
ment characteristics of 9l c which are entirely similar to those imposed on 
91* by condition (2). 

(iv) Finally, the elements of 9l B and 91° satisfy analogous Lipschitz 
conditions by condition (4). 

Therefore, the variables (e e ,q c ) and (E c ,q) of 9l c are analogous to 
the variables (e r , i r ) and (E*J*) of 91*. 

Remembering that E c is linearly related to E, we conclude that the 
voltages and charges of the capacitors at time t are uniquely determined 
by the values, at the same time t, of the voltage sources E (0 and q (/), 
and that the mapping 

e B (t) = f, c (E(0,q(0,0 (10) 

q f (/) = f„ e (E(*),q(0,0 (ID 

satisfies Lipschitz conditions. 

Since #, in any fundamental cut set is equal to the sum of the capacitor 
charges 

'h£l = jnt) 

at 

where ./, c (0 is the contribution of the current sources to the ith cut set. 
As J c is a linear combination of i r , U , and J it follows that 
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|q(0 = f 3 (ir(0,u(0,J(0) (12) 

where f q is linear and does not depend explicitly on time. 

Let us now consider the inductors of the network 91. Given any set of 
resistor voltages e r and capacitor voltages e c such that Kirchoff's voltage 
law is satisfied in each loop formed by capacitors, resistors and voltage 
sources, let us replace the capacitor and resistor branches of 91 with volt- 
age sources equal to the corresponding voltages e c , e r . The network now 
consists of sources and inductors only. Let us again use the algorithm of 
Section III to change the configuration of the sources and separate the 
network into its separable parts. The resulting network is denoted by 
9l 7 ' and its sources by E ,J '. As in the case of the resistive network, we 
are going to use the result previously proved for 91" to deduce a similar 
result for 91 . 

91 L consists of the three sets of subnetworks which were described in 
connection with 91*. Consider the third set of subnetworks of 9l L , which 
consists of single inductors in parallel with a current source. Except for 
the trivial case where 91 consists only of a single inductor in parallel with 
a current source, this set is empty, for otherwise Tfl(E)u)' would be sepa- 
rable. Condition (3) implies that each subnetwork of the first set has a 
tree, say t l , whose tree branches are monotonically increasing and cur- 
rent-controlled and whose links are monotonically increasing and flux- 
controlled. For each subnetwork of the first set, with each fundamental 
loop of t l defined by an inductive link, we assign a variable <p t equal to 
the sum of the fluxes of all the inductors of that loop. For each subnet- 
work of the second set we assign a <pi equal to the flux of the inductor. 
«£ will denote the vector whose components are the ^,'s. 

The analogy between 91 ^ and 91 ff is established in four steps: 

(t) For 91 w the Kirchhoff voltage law applied to the tth loop associ- 
ated with a resistive link reads 

k 

where e*< is the voltage across the feth branch of the tth loop. The e ki 
are components of e r . For 91 L we have by definition of <pi , 

<Pi = J! <Pki 
k 

where ipui is the flux in the fcth branch of the ith loop. The <p ki are com- 
ponents of <q l . 

(it) For both 91 R and 91 L the Kirchhoff current law holds. 

(m) Condition (3) imposes requirements on the topology and ele- 
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ment characteristics of 91 '' which are entirely similar to those imposed 
on 91* by condition (2). 

(iv) Finally, the elements of 91 " and 91 *" satisfy analogous Lipschitz 
conditions by condition (4). 

Therefore, the variables (tQ L , i L ) and (<qJ l ) of 91*" are analogous to 
the variables (e r , i r ) and (E J 1 ), 

As J L is linearly related to J, we conclude that fluxes and currents of 
the inductors at time / are uniquely determined by the values, at the 
same time t, of («e(OJ(0)> anc ' tnat the mappings 

v*(0 - **., MO J (0,0 (13) 

U(0 -fu(v(0J(0^) (14) 

satisfy Lipschitz conditions. 

Since y>, in each loop is equal to the sum of the fluxes in the loop, it 
follows from the Kirchhoff voltage law that 

dt 

where E* {t) is the contribution of the voltage sources in the ith loop. 
As E L is a linear combination of e c , e r , and E, it follows that 

^s?(0 = f,(e o (0,e r (0,E(0) (15) 

where f v is linear and does not depend explicitly on time. 

Any solution of the network requires that (8), (9), (10), (11), (12), 
(13), (14) and (15) be satisfied simultaneously. It is shown in the follow- 
ing that these equations determine a unique solution. 

In (12) and (15) substitute values of e c , e r , i r and \ L from (10), (11), 
(9) and (14). The results are 

j t q - Ufi r (i c (E A ,t),E,i iL (^J,t), J,t),t iL (v, J,0, J] (16) 

!«<>= f,[f ( ,(E,q,/,),f, r (f, c (E,q,/),E,f 1 J^,J,n,J,0,E]. (17) 

Since the right-hand sides of (16) and (17) are compositions of func- 
tions satisfying Lipschitz conditions, these equations may be rewritten 
as 

|q(0 = F a (E(0,q(0^(OJ(0,0 (18) 
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i v (0 - f,(e(0a(0*(0J(0i0 < 19 > 

where F g and F v satisfy Lipschitz conditions in q and <q. Therefore, for 
any E(-) and J(-) that are regulated functions of time 13 and for any 
initial values of cp and q, the differential equations (18) and (19) deter- 
mine uniquely <e( • ) and q(-), and the solutions are continuous functions 
of time. 18 In terms of E(-), ](•),<&(■) and q(-), equations (8), (9), (10), 
(11), (13) and (14) determine uniquely the currents and voltages of the 
resistive branches, the voltages and charges of the capacitive branches, 
and the fluxes and currents of the inductive branches. Therefore the 
network 91 is determinate. (Incidentally, the proof shows that the state 
of the network may be represented by (q, <£»).) 

It is worth indicating an immediate consequence of (18) and (19) and 
the other circuit relations. 

Corollary: If the conditions of Theorem IV are satisfied, E and J are 
continuous functions of time, and all elements depend continuously on 
time, then e c , q, , e r , l , <g L , U are continuous functions of time; in 
other words, jump phenomena 1 are excluded. 

Corollary: Let the network 91 consist of independent sources, nonlinear 
(possibly time-dependent) nionotonically increasing one-to-one re- 
sistors, capacitors and inductors. If the characteristics of all elements 
Lipschitz conditions as described in condition (4) of Theorem IV, the 
network 91 is determinate. 

Corollary: If branches of a network 91 consist of: (a) voltage sources, 
current sources; (b) one-to-one nionotonically increasing resistors, 
capacitors, and inductors whose characteristics satisfy condition (4) of 
Theorem IV; (c) one-to-one two-poles of the types described by the 
theorems of Section II and which satisfy conditions (a), (6), (c) and 
(d) of these theorems: then network 91 is determinate. 

Given a physical circuit or device, it may happen that a particular 
model of the circuit does not satisfy the conditions of Theorem IV. For 
example, this model 9l' might be such that S*fSl\nc includes a parallel 
connection of two charge-controlled capacitors, Di(q) and D 2 (q), with 
only Di monotonically increasing. Under these conditions, it may happen 
that the current through the parallel combination does not determine 
uniquely the voltage across it. If, however, the model is changed (call it 
91") and a resistor (or inductor) is inserted in series with D 2 , then 
S*9l"(*)c now includes an open branch D 2 , condition (1) of Theorem 
IV is no longer violated, and 91" is determinate. Obviously, this idea 
may be used in the case of inductors and resistors. 
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Finally, let us conclude this section by a discussion which draws at- 
tention to some consequences of the conditions of Theorem IV. Some of 
the properties considered hero will be used in the next section for writing 
in detail the network equations. 

Let 91 be a network which satisfies the conditions of Theorem IV. De- 
note by ?yi_i , g the sets of resistors which are current-controlled (but 
not voltage-controlled), one-to-one, and voltage-controlled (and not cur- 
rent-controlled), respectively. Similarly, denote respectively by d,di-i , c 
the charge-controlled (and not voltage-controlled), one-to-one, and volt- 
age-controlled capacitors, and by /,rj_i , T the current-controlled (and 
not flux-controlled), one-to-one, and flux-controlled inductors. 

Let us carry out the following operations: 

(a) choose a forest of 91 B 

(b) choose a forest of 9l( B ),/ 

(c) choose a forest of 9l( B( /),i 1 _ 1 

(d) choose a forest of 9l(Md,_ 1 ) r 

(e) choose a forest of 9l(K<w 1 _ 1 r)r i _ 1 

(f) choose a forest of 9l<i?<M,_ 1 rr,_ 1 )/ 

(g) choose a forest of 9l<AWi_ 1 rr,_ii)ri_i . 

Since the conditions of Theorem IV are satisfied by 31, it follows that 
the union of these forests forms a tree of 91 which we denote by t. The con- 
struction of this tree is an extension of Bryant's procedure. 

This can be proved in the following way: From the conditions of 
Theorem IV it follows that the union of the forests chosen by (b) and 
(c) [by (d) and (e), by (f ) and (g)] are forests of JJl^c , [9l<*o* , and 
SliBCML, respectively]. Let us add the network's resistors to 91(bc«)l- 
This is done by splitting nodes and adding the new branches between 
them. Consider a node which was split, say, to three nodes and a resistor 
subnetwork connected between these nodes, It is clear that the subtree of 
this resistive subnetwork completes the forest of 31 {B cr)l for a forest of 
the yi ( EC)RL • We can use the same argument to show that by adding the 
capacitors and voltage sources we get a forest of the network which in- 
cludes all branches but the current sources. However, the current sources 
do not form any cut set and therefore are links of this forest. Thus t, 
the union of these forests, is a tree of 91. 

From the construction of the tree, the conditions of Theorem IV, and 
the above discussion it follows that r contains all current and charge-con- 
trolled elements which are not one-to-one, and all voltage and flux-controlled 
elements which are not one-to-one are links of this tree. 

Consider the fundamental cut set of 91 defined by an element of set d, 
a charge-controlled capacitor whose characteristic is not monotonically 
increasing. By assumption, this capacitor is an open branch of S**3l {B )c ; 
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this cut set may not contain capacitors or voltage sources and therefore 
consists solely of resistors, inductors and current sources. Similar proper- 
ties exist for fundamental loops defined by the various links. One can ex- 
hibit these properties by making a table in which links and tree branches 
are partitioned according to the types of their elements and properties 
of their characteristics; for each link and branch the table specifies the 
type of elements that are allowed to be in the corresponding loop or cut 
set. As the table is complicated, it is omitted and only some of the more 
interesting properties are listed below. Here we are going to make use of 
the rank-order of the elements, ECRLJ, defined in Section II: 

(i) Tree branches with characteristics which are not monotonically 
increasing are the highest ranked elements in their own fundamental cut 
set. Thus, for example, a charged-controlled nonmonotonically increas- 
ing capacitor has a fundamental cut set which may include links which 
are resistors, inductors and current sources but no other capacitors. 

(ii) Links with characteristics which are not monotonically increasing 
are the lowest ranked elements in their own fundamental loop. Thus, a 
fundamental loop defined by a nonmonotonically increasing resistor 
may have only capacitors or voltage sources in its tree branches. 

VI. EQUATIONS FOR RLC NETWORKS 

The purpose of this section is to write explicitly the equations of a 
nonlinear RLC circuit of the type considered in (he previous section. 
Another purpose is to exhibit the similarities and differences between 
the equations that describe linear networks and those that describe the 
nonlinear networks under consideration. 

To simplify the exposition consider the resistive network of Fig. 14. 
Call t\ the tree formed by the branches 1,2,3, and the voltage source E. 
If the network were linear, the fundamental cut set equations would 
read 




Fig. 14 — Resistive network. 
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Q#\ + Q\ifi\ + e 2 — e 3 — E) = J i 

04 (ci + e 2 - e 3 - E) + 2 e 2 - £ 5 (-e 2 + e a ) = /i (20) 

-04 (ei + e 2 - e 3 — E) -f 03^3 + B (-e 2 + e 3 ) = 

where e,- is the voltage across the tth branch and 0, is the conductance 
of this branch. In the well known matrix form, the equations become 
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(21) 



or, more generally, 

At(R),rGr& T(RE).R -p 



+ A 



T(R).UJ) 



J = o 



(22) 



where e, E and J are column vectors whose components are the tree 

V 

the voltage sources, [E], and the current sources, 



branch voltages, 



[./], of the network; G H is the branch admittance matrix. The A's are 
appropriate submatrices of the fundamental cut set matrix Q. The first 
subscript of A denotes the rows and the second subscript denotes the 
columns of Q whose intersection forms the submatrix. Thus A t(RB ),l(r) 
is a submatrix formed by the intersection of rows corresponding to 
resistive and voltage source tree branches and columns corresponding 
to resistive links; Ar(«).« is formed by the intersection of rows corre- 
sponding to resistive tree branches and columns corresponding to re- 
sistive branches. A r («>,j is defined similarly. The prime over a matrix 
indicates transposition. Now, let the resistors become monotonically in- 
creasing one-to-one nonlinear resistors. Without loss of generality we 
can assume these new resistors to be time invariant. Let 0i(- ), 2 (-)> 
03 (Oi 04(') and 0s(-) be their characteristics. 
The cut set equations are : 

0i (d) + 04 (ei + eo - e 3 - E) = Ji 

fc(d + c-2 - c 3 - E) + 2 (e 2 ) - 5 (e 3 - e 2 ) = /, (23) 

-04 (ci + Bi — ez — E) + 05 (ea - e 2 ) + 03 (e 3 ) = 

where, for example, g\(t ; \) is now the value of the function g x evaluated 
at e\ . 

The similarity between (20) and (23) suggests a shorthand notation 
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for writing the equations of nonlinear networks. By the product A*x 
(where A is a diagonal matrix whose elements are functions a, ( • ) and 
x is a column vector whose components are xi , x 2 , • • • , x n ), we denote 
the column vector whose ith. component is a,(rc,), that is, the «'th di- 
agonal element of A evaluated at the ith component of x. With this 
symbolic notation the equations of the network of Fig. 14 can be written 
(for the nonlinear case) in a form analogous to (22). 



'T(ff) 



,rGr* (A T{RB),R J? ) 



+ &T(R),L(J)J — 



where G H is the diagonal matrix whose elements are the characteristics 
Q\ > Qi i • ' ' > 06 and the * operation must be interpreted as indicated 
above. G R will be referred to the branch characteristic matrix. With this 
symbolic notation, cut set matrices, loop matrices and branch resistance 
matrices may be used to writing equations of nonlinear networks in the 
same way as for linear networks. 

Let us now assume that the elements of the tree n of the network of 
Fig. 14 are monotonically increasing current-controlled but not voltage- 
controlled and the links are monotonically increasing voltage-controlled 
but not current-controlled. Since the tree branches are not voltage- 
controlled, the equations cannot be written in the form of (22). Let 
h(')i ft(') and f 3 (- ) represent the characteristics of the tree branches 
and ii , i* and i 3 be the currents of the corresponding tree branches. In 
terms of the tree branch voltages and currents the cut set equations 
become: 

i\ + <}\ (ei + e 2 — e 3 — E) = J x 

ii + Oa (ei 4- ea + e 8 - E) - fo(e 3 - e%) = J\ (24) 

i-.\ - 0i(ei + e 2 — e 3 — E) + g & (e 2 - e s ) = 0. 

The other set of equations is 

^i = n(ii) 

e 2 = fate) (25) 

c 3 = fsiis) 

or symbolically 

ir(ff) + &T(R),L(R)Gl(R)* I A T(RE).UR)\ p ) 4" At(R),L(J)J = 0. (26) 

er(R) == Rt(R)*1t(r) (27) 

I 
where i T (R) , b T (r) , E and J are the tree currents and voltages, and vol- 

/ 
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tage and current sources respectively. G and R are the link and branch 
characteristic matrices, and in our example 



G L (R) — 
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The A's are appropriate sub-matrices of the fundamental cut set matrix 
Q. The first subscript denotes the rows and the second subscript denotes 
the columns of Q whose intersection forms the submatrix. Thus 
A T (s«),L(B) is a submatrix formed from the intersection of rows cor- 
responding to resistive and voltage source tree branches and columns 
corresponding to resistive links. &TUt),um and &t(r),hj) are defined 
similarly. A comparison of (20), (27) and (22) shows that in the case of 
current-controlled tree branches and voltage-controlled links which are 
not one-to-one, we need both ir<«) and er(«) for a straightforward writ- 
ing of the cut set equations and the branch characteristic equations. 
Either ir(R) orer<«) can be eliminated from the equations. The resulting 
equations are: 



l7 - (H) T ^r(R),UK)G i,(li) : 



Or 



T(RE),L(R) 



R 



T(R)*lT(R) \\ 

E \) 



(28) 



+ A 



T{R),L(J) 



J = 



&T(K) "T" Rt(R)* \ A !•(«), MR)Gl(R)* I A T(RE),L(R) 



er<«) 

E 



]) 



(29; 



+ R 



T(«) ! 



(A 



T(R).HJ) 



J = 0. 



Fundamental loop equations can be written in a similar way using both 
the voltages and currents of the links i/, ( «) and e /.</,•> . The equations are 



■L(R) 



+ Il(R),T(.R)Rt(R)* [ 1 URJ).T(R) 



['?']) 



+ 1 



UR).T(E) 



E = 



iz-(ff) — Gi 



,*e 



Z.(W)*C/,(«) 



where the l's are appropriate subniatrixes of the fundamental tie set 
matrix B. Similarly to (26) and (27), cither e L ,K) or i/. (W) can be elimi- 
nated. 

We now write the equations for a general RLC network (which satis- 
fies the requirements of Theorem IV) by performing the following steps f 

t Other systems of variables are possible. For example, one can choose charges 
and fluxes as above and voltages of resistive links whose loop does not consist of 
capacitors and voltage sources only. 
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(i) A tree is chosen as explained in Section V. 

(it) Variables are chosen. We choose here the charges on the capaci- 
tive tree branches, q D , the currents of the resistive tree branches whose 
fundamental cut set does not consist of inductors and current sources 
only, i R , and the fluxes of the inductive links, <pr . 

The equations make use of the following characteristic branch ma- 
trixes: C and D are diagonal matrixes whose elements are the charac- 
teristics of voltage-controlled and charged-controlled capacitors, re- 
spectively; G and R, those of voltage-controlled and current-controlled 
resistors, respectively; L and r, those of current-controlled and flux- 
controlled inductors, respectively. Without loss of generality we can 
assume that the elements are time-invariant. The equations are 



dt 



« _L A C *( A' nDr<c)*q D T\ 

q D "f &T(C),L(C)l'L(C)*i A T(CB),t(C) jj, J 



-f- Aj-(C) , L(R)G h(lt)* ( A T(CKE),L 



(K) 



D T (D*qD 

Rt(r)*1r 

E 



(30) 



+ AjT(C),L(i)rL(i)*<ffr + ^T(C),L(J)J — 



1« T ^T(R) , L(R)G L(U)* I A T(CRE),L(R) 



+ A 



T(R).L(L)1 L(L) 



Dt(c)*<Id 

Rt(R)*1r 

E 



r««*(5r + A 



(31) 



T(R),L(J) 



J = 



P * /V rrz.(L)*«grl\ 

RO K nRi)* \ l L(LJ),T{R X ) j N 



/.(t).T(i 



(32) 



"t - iu.1.) ,t(c)Dt(c)*Qd ~r ll(l),t(r)Rt(r)*1r 

"T 1l(L),7'(B)E = 

where #1 is the set of resistive tree branches whose fundamental cut set 
contains inductive links and current sources only; and R is the set which 
contains all other resistive branches. 

The terms in the brackets in (30) and (32) are equal to our state 
variables q and «p of Section V. One can write the equations in terms of 
these variables : the relations between q D and q and <£r and <j> are given by 



NONLINEAR RLC NETWORKS 197 

= q 



<1D T ^T(C),L(C)^UC)* I ^* T(CB).L(.C)\ p, I 

.„ 4- 1 T ±(\' \ T L{L)*iQT~\\ 

tffr ~T lt(L),T(L)*T(£)* I I L(LJ),T(L)\ •, I 



- <0- 



In summary, the equations of the /?LC nonlinear network are written 
in a way which is a generalization of the methods used in linear networks. 
However, great care must be taken of the fact that some characteristics 
are representable by functions w hich do not have inverses. This section 
indicated a method for tackling the problem. In this section, the equa- 
tions are written in terms of three sets of variables: q D , the charges on 
the capacitive tree branches; yg L , the fluxes in the inductive links and 
iff , the currents in the resistive tree branches whose fundamental cut 
sets do not consist of only inductors and current sources. It is interest- 
ing to note that (except for the trivial case where 91 consists of a single 
capacitor in parallel with a voltage source or a single inductor in parallel 
with a current source) the dimension of the state vector (q,<£) used 
above is the same as in the linear case: [number of independent initial 
conditions] = [number of reactive elements] — [number of independent 
capacitor-only tie sets] — [number of independent inductor-only cut 

j. 1 7,8,9 
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